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ABSTRACT 

Drawing on recent estimates of the power of jets from X-ray binary systems as a func- 
tion of X-ray luminosity, combined with improved estimates of the relevant Log(N)- 
Log(Lx) luminosity functions, we calculate the total energy input to the interstellar 
medium (ISM) from these objects. The input of kinetic energy to the ISM via jets is 
dominated by those of the black hole systems, in contrast to the radiative input, which 
is dominated by accreting neutron stars. Summing the energy input from black hole 
jets Lj in the Milky Way, we find that it is likely to correspond to > 1% of LsNe, the 
time-averaged kinetic luminosity of supernovae, and > 5% of Lqr, the cosmic ray lu- 
minosity. Given uncertainties in jet power estimates, significantly larger contributions 
are possible. Furthermore, in elliptical galaxies with comparable distributions of low 
mass X-ray binaries, but far fewer supernovae, the ratio Lj/LsNe is likely to be larger 
by a factor of ~ 5. We conclude that jets from X-ray binaries may be an important, 
distributed, source of kinetic energy to the ISM in the form of relativistic shocks, and 
as a result are likely to be a major source of cosmic rays. 
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The interstellar medium (ISM) in the Milky Way and other 
galaxies receives a continuous input of energy in the form of 
both bulk motions and radiation. The result of this energy 
input is the heating of the ISM, the generation of interstellar 
turbulence and the production of high energy cosmic rays. 
The most violent processes are those associated with stel- 
lar evolution: stellar winds, protostellar jets and supernovae 
(SNe). More specifically, cosmic ray production is only likely 
to be associated with high velocity (i.e. J; 0.1c) shocks, and 
so their most likely source has long been considered to be 
the SNe. See e.g. Chevalier (1977), McCray & Snow (1979) 
and Elmegreen & Scalo (2004) for reviews of these processes. 

However, in recent years it has become clear that there 
is at least one other potential source of powerful relativistic 
shocks in the ISM, namely the jets resulting from accretion 
onto stellar mass black holes and neutron stars in X-ray 
binary systems (Mirabel & Rodriguez 1999; Fender 2005). 
These jets seem to be powerful and ubiquitous in the ma- 
jority of such systems which exist in 'hard' X-ray spectral 
states and/or are undergoing X-ray outbursts (Fender 2001; 
Fender, Belloni & Gallo 2004). Heinz & Sunyaev (2002) es- 
timated the total power input from a handful of the most 



powerful jets in our galaxy and concluded that they might 
contribute 1-10% of the cosmic ray luminosity of the galaxy, 
perhaps with a specific spectral signature. However, recent 
advances in our estimates of the scaling of jet power with X- 
ray luminosity (Fender, Gallo & Jonker 2003; Fender, Belloni 
& Gallo 2004), combined with improvements in our knowl- 
edge of the distribution of X-ray sources as a function of 
luminosity (e.g. Grimm, Gilfanov & Sunyaev 2002,2003; Gil- 
fanov 2004) have allowed us to recalculate more accurately 
the sum of the power input from such jets. In this paper we 
present the results of these calculations. 



2 LUMINOSITY FUNCTIONS 

2.1 Differential luminosity function for LMXBs 

Grimm et al. (2003) presented the X-ray luminosity function 
(XLF) for LMXBs in the Milky Way in the form of a cut- 
off power law. Since then, Gilfanov (2004) has produced an 
average differential XLF for LMXBs in all galaxies. This 
function is in the form of a power law with two breaks: 
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Figure 1. Estimated total numbers of X-ray binaries in the Milky 
Way as a function of X-ray luminosity. At the highest luminosi- 
ties, above 10 35 erg s _1 , the low-mass X-ray binaries (LMXBs) 
dominate the source population; although at even lower luminosi- 
ties numbers of high-mass X-ray binaries (HMXBs) may come 
to dominate, the numbers are not well measured at those levels. 
Since the process of disc-jet coupling does not seem to be strongly 
affected by the nature of the companion star (see e.g. discussion 
in Fender 2005) this figure illustrates that the total X-ray binary 
jet power into the interstellar medium will be dominated by the 
LMXBs. 
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where Lx,38 = Lx/W 3 
are related by 



erg s 1 and normalizations ^1,2,3 



Ki = Ki (L 6 ,!/L M ) a2 
K z = K2(L b , 2 /L cut ) a3 

The best fit to the overall normalisation is given as 
Ki = 440.4 ± 25.9 per 10 n M Q (2) 



In the following we shall use this as the template LMXB 
XLF for all galaxies. Based on Liu et al. (2000, 2001) we 
estimate the black hole and neutron star fractions amongst 
the LMXBs to be f BH = 0.2, f NS = 0.8. This luminosity 
function, normalised to the Milky Way using a total stellar 
mass of 4.5 x 10 10 M Q (Gilfanov 2004) is plotted in Fig 1. 



2.2 Differential luminosity function for HMXBs 

For HMXBs, we use the differential form of the XLF from 
Grimm et al. (2003), which is simply a steep power law: 
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(3) 



where for the Milky Way K ~ 0.7. Based on Liu et al. 
(2000, 2001) we estimate the black hole and neutron star 
fractions amongst the HMXBs to be /bh = 0.05, /jvs = 
0.95. This function is plotted in Fig 1 alongside the LMXB 
XLF. 



2.3 Lj(Lx) for X-ray binary jets 

2.3.1 Steady ('low /hard' state) jets 

Corbel et al. (2000; 2003) and Gallo, Fender & Pooley 
(2003) have found an apparently universal correlation be- 
tween radio and X-ray luminosities for black hole candi- 
date (BHC) binaries in the 'low/hard' and 'quiescent' X-ray 
states (which seem to correspond to accretion rates below 
about 1-10% of the Eddington limit). The relation has the 
form 



radi 



ocix 



(4) 



where b ~ 0.7. The relation between jet and radio lu- 
minosities may also be considered to have a power-law form 



oc Lj 



(5) 



Several models of steady, conical jets predict c ~ 1.4 
(e.g. Blandford & Konigl 1979; Heinz & Sunyaev 2003), 
which leads to 



Lj — Asteady.BHcLx 

in Eddington units, or 



JJ.38 
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steady, BHC 



(Lx,3 



M 



(6) 



(7) 



indicating that the fractional jet power increases rapidly 
as the X-ray luminosity decreases (Fender, Gallo & Jonker 
2003; Fender, Belloni & Gallo 2004). 

The value of the normalisation A tea dy,BHC i s uncertain. 
Fender, Gallo & Jonker (2003) estimated (in their opinion 
conservatively) that ^steady, bhc > 6 x 10~ 3 . While it is 
certainly far from the consensus, larger values of A s t e ady,BHC 
may be more realistic. Malzac, Merloni & Fabian (2004) have 
suggested that even at Lx ~ 10~ 3 L E dd, the jet may be an 
order of magnitude more powerful than the X-ray emission, 
which corresponds to A tea dy,BHC ~ 0-3- Alternatively, we 
can set the value of A stea dy,BHC so that it corresponds to 
the transition from the low/hard to high/soft X-ray states 
(see e.g. McClintock & Remillard 2005 for a discussion), in 
which case Abhc ~ 0.1. In all the following discussions we 
shall consider 0.006 < Abhc < 0.3 as covering the full range 
of likely values. 

We further assume that A stea d y ,NS = 0.1A stea dy,BHC, 
based on the lower 'radio loudness' of NS XRBs (Fender & 
Kuulkers 2001; Migliari et al. 2003; Fender, Gallo & Jonker 
2003; Muno et al. 2004; Migliari & Fender in prep). It should 
be stressed that the Lradio oc L^' 7 relation has not been es- 
tablished for neutron star X-ray binaries, and is only as- 
sumed in Fender, Gallo & Jonker (2003) by analogy with 
the BHCs. However, since the BHCs are so much more radio 
loud then this is not significant for the total energy budget 
calculated here. 



2.3.2 Lj(Lx) for transient jets 

A fit to the power in transient optically thin ejection events 
from black hole X-ray binaries as a function of X-ray lu- 
minosity (Fender, Belloni & Gallo 2004), simplified to one 
significant figure, gives 



Lj — A 



T 0.5 

trans, BHC-^X 



(8) 



© 0000 RAS, MNRAS 000, 000-000 



3 



in Eddington units, where A tr ans,BHC ~ 0.4. The same 
L x ' 5 dependence (see Fender, Belloni & Gallo 2004 for 
caveats) as for the steady jets means the equation has the 
same form with simply a different normalisation: 



JJ.38 



Kr.u Bl)c ' L s ,3s)°' J ( -^j— 



(9) 



It should be stressed that the normalisation is also 
rather uncertain, being based on a small sample of sources 
for which the energy estimate is itself intrinsically uncertain 
and different to that used for the steady jet sources (see dis- 
cussion in Fender, Belloni & Gallo 2004). Nevertheless, we 
consider it likely that this estimate is accurate to within one 
order of magnitude, and that if it is not then it is more likely 
to be an underestimate than an overestimate. Based upon 
this we will consider in this paper the range of likely values 

tO be 0.04 < Arans.BHC < 4.0. 

For the transient sources we further need to consider 
periods in 'soft' X-ray states when the radio emission, and 
hence presumably the jet, are suppressed or 'quenched' 
(Fender et al. 1999; Corbel et al. 2001, 2004; Fender, Belloni 
& Gallo 2004; Fender 2005). The physical processes responsi- 
ble for the suppression of the jet are not clear, but they seem 
to reduce the radio emission by a factor > 30 - as a result we 
shall consider that when in such states the jet power drops 
to zero. The duty cycle in such states varies from outburst 
to outburst and it is hard to estimate a reasonable mean 
value. Based upon outbursts of the black hole LMXBs GX 
339-4, XTE J1550-564 and XTE J1859+226 it seems a value 
of 60-80% for the duty cycle in 'soft' jet-suppressed states 
is reasonable (J. Homan, private communication; see also 
Fender, Belloni & Gallo 2004; Homan & Belloni 2005). As a 
result, we multiply the integrated power from transients by 
a factor 0.3. While this number is rather uncertain, we note 
that it is still the uncertainty in the value of the normalisa- 
tion Arans which dominates the integrated power input. 

In the absence of sufficient good data on neutron star 
transient outburst with which to make a fit to A tr ans,NS, we 
simply follow the approach adopted for the steady jets, i.e. 



itrans,NS 



O.lAr 



;,BHC- 



3 TOTAL POWER 

As noted above, -Pj(Lx) is probably a function both of X-ray 
state and of the nature of the accreting compact star (neu- 
tron star or black hole). Convolved with both HMXB and 
LMXB luminosity functions, this presents us with eight pop- 
ulations of jet producing objects to consider (LMXB/HMXB 
x steady /transient x NS/BHC). 

The integrated total jet power for any one of these 
classes is given by 



J ,1, total 



1X,: 



D dN AT 



(10) 



where Pj(Lx) is the jet power at a given X-ray luminos- 
ity Lx, and NxdLx is the differential luminosity function for 
the class of object. The term / is a correction factor which 
can be separated into two parts: 



/ — fpopfs 



(11) 



The first term, f pop is simply the fraction of the pop- 
ulation of that class of objects, in the sense that f pop ,NS + 
f pop ,BHC = I- The second term, f spe c compensates for the 
fact that the X-ray luminosity functions given above, from 
Grimm et al. (2003) and Gilfanov (2004) consider the X-ray 
luminosity in the range < 10 keV, whereas the expressions 
for the jet power were arrived at considering the total X-ray 
power. This term is therefore in some sense a bolometric 
correction. The values of f ap ec for black holes are estimated 
in Portegies Zwart, Dewi & Maccarone (2004) for the black 
hole candidates. For the neutron star sources we followed 
essentially the same method. Estimated values for f pop and 
/spec are presented in table 1. 

It is already clear from Fig 1 that HMXBs are unlikely 
to contribute significantly to the total power output, kinetic 
or radiative, of X-ray binaries (with the notable exceptions 
of SS 433 and Cyg X-3 which are powerful jet sources; see 
section 4). Simple calculations bear this out and from this 
point on we do not consider HMXBs further. 

Radiative and kinetic output of LMXBs 

In Fig 2 we present the integrated X-ray and jet (kinetic) 
luminosities of LMXBs as a function of source X-ray lumi- 
nosity, scaled for the Milky Way, for three sets of values of 
jet power normalisations A stoa( jy and A tr ans. The shape of 
the functions result from the multiplication of the LMXB 
XLF (equation (1)) and the switch from steady to transient 
jets, which we fix to occur at 1.9 x 10 37 for the black holes 
and 6 x 10 36 for the neutron stars, e.g. around 20% of the 
Eddington limit for each class of object. Some features are 
common between the figures, and indeed independent of the 
selected values of the jet normalisation, namely: 

• The radiative (X-ray) output is dominated by the more 
numerous NS LMXBs 

• The kinetic (jet) output is dominated by the BHC 
LMXBs, which produce more jet power at a given Lx 

The integrated bolometric radiative (X-ray) output of 
the two classes of objects is L x ,boi,N,s — 6.4 x 10 39 erg s _1 
and L x ,bol,BHC = 1.7 x 10 39 erg s -1 . 

Moving on to individual figures, Fig 2(a) ('case A') con- 
siders the 'minimum' or conservative estimates of the jet nor- 
malisation. In this scenario both classes of object, NS and 
BHC, produce more integrated radiative than kinetic ouput. 
Even with these low values for the normalisation, a transi- 
tion from 'X-ray dominated' to 'Jet dominated' is apparent 
for the BHCs (see Fender, Gallo & Jonker 2003; Malzac, 
Merloni & Fabian 2004 for a fuller discussions). The inte- 
grated BHC jet luminosity is Lj = 1.5 x 10 38 erg s _1 (22% 
from transient jets). 

Fig 2(b) ('case B') considers the situation where the 
steady jets reach power equipartition with the X-rays at 
the state transition, and dominate at lower X-ray luminosi- 
ties (see e.g. Livio, Pringle & King 2003 for support for this 
point of view). The transient jets have a normalisation corre- 
sponding to the best fit from Fender, Belloni & Gallo (2004) . 
In this scenario the peak of the kinetic power output from 
the BHC LMXBs is comparable to the peak of the radia- 
tive power output from the NS LMXBs. The integrated jet 
luminosity is Lj = 2.0 x 10 39 erg s _1 (19% from transient 
jets). 
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Table 1. Estimated values for the population f p op and spectral / s pec corrections to be made to the integrals for the different classes of 
X-ray binary, as well as the jet power normalisation A. See text for more details. 
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Table 2. Supernova rates in SNu for different galaxy types 
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Figure 2. The summed power output of all the Milky Way X-ray 
binaries, in radiation (X-rays) and kinetic power (jets), for three 
different sets of jet power normalisations, as a function of soft X- 
ray luminosity. The neutron stars dominate the radiative output 
but the kinetic power output is dominated by the black holes. 



Galaxy type 
E-S0 
SOa-Sb 
Sbc-Sd 



la 

0.18 ± 0.06 
0.18 ± 0.07 
0.21 ± 0.08 



Ib/c and II 

< 0.03 
0.53 ±0.19 
1.00 ±0.35 



Total 
0.18 ±0.06 
0.72 ±0.21 
1.21 ±0.37 



Fig 2(c) ('case C) is perhaps the most 'extreme' case, in 
which the kinetic output from the BHC LMXBs dominates 
the radiative output at essentially all Lx- The integrated jet 
luminosity is Lj = 1.0 x 10 40 erg s _1 (26% from transient 
jets). 



4 JETS VS. SUPERNOVAE 

A 'standard' type-I or type-II supernova explosion is be- 
lieved to deposit around 10 51 erg of kinetic energy into the 
surrounding interstellar medium (e.g. Chevalier 1977; Korpi 
ct al. 1999). For an estimated rate of one such event per 100 
years in the Milky Way, the time-averaged kinetic power is 
LsNe = 3 x 10 41 erg s _1 . Furthermore, Heinz & Sunyaev 
(2002) estimate the cosmic ray luminosity of the Milky Way 



as being approximately 10% of this, Lcr = 4 x 10 erg s~ . 

Comparing these numbers with the three cases of jet 
power normalisation considered above, we can see that 

10~ 3 <, Lj/LsNe ^ 0.03 

10~ 2 <, Lj/Lcr <, 0.3 

So it is clear that there may be enough power in the 
jets from X-ray binaries to supply a sizeable fraction of the 
cosmic ray luminosity of the Milky Way. There are many 
caveats to this, of course, but this number is of considerable 
interest. If we consider 'case B' outlined above to be the 
most likely, we find that the summed kinetic luminosity of 
X-ray binary jets is <; 5% of Lcr. 



4.1 Early- type galaxies 

So far we have only considered the situation for the Milky 
Way. However, different types of galaxies have different ra- 
tios of X-ray binaries to supernovae, which will affect the 
ratios calculated above. Estimated supernova rates for dif- 
ferent types of galaxies are given in table 2 (from Cappellaro, 
Evans & Tuatto 1999). 

Note that the total SNe rate in elliptical galaxies is 
about a factor of six less than in late-type spirals such as 
the Milky Way (due to the complete absence of SNe types 
Ib/c and II - see e.g. Cappellaro et al. 1999; van den Bergh, 
Li & Fillipenko 2002). 

What do we expect of the X-ray binary, jet-producing, 
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population ? Gilfanov (2004) has shown that the popula- 
tion of LMXBs scales linearly with the stellar mass of the 
galaxy, whereas for early-type galaxies the HMXB popula- 
tion is expected to be essentially zero. We should note that 
at least two of the most powerful jet sources, Cyg X-3 and 
SS 433, appear to be HMXBs (but in fact were not included 
in the preceding analysis as all HMXB contributions were 
disregarded). Nevertheless, GRS 1915+105 is a LMXB as 
are the seven most luminous neutron star XRB jet sources 
(the Z sources plus Cir X-l), as are all the BHC transients. 
Based on this we may crudely estimate that the jet power 
from 'outbursting' systems, per unit mass, is reduced by 
about a factor of no more than two in ellipticals compared 
to late-type spirals. 

As an example we may take the case of NGC 4472, 
a large early-type galaxy. Gilfanov (2004) gives its stellar 
mass (the component which is correlated with the number 
of LMXBs) as being 3.5 times that of the Milky Way. Con- 
sidering the scaling of LMXB numbers with stellar mass, as 
well as the possible 50% reduction in the brightest sources 
due to the lack of 'odd' HMXBs such as SS 433 and Cyg 
X-3, we would expect the summed jet power in NGC 4472 
to be about twice that in the Milky Way. The total number 
of supernovae should also scale by mass, but is reduced by 
a factor ~ 6 due to the lack of any type Ib/c or type II 
events, and so would be around half of that in the Milky 
Way. Therefore the ratio of Lj/Lsnb would be ~ 4 times 
greater in this system, which may be considered to by typical 
for early- type ellipticals. In this situation, 'case C outlined 
above would correspond to a total jet power input which was 
*J 15% of the power output of SNe. Even in the more con- 
servative 'case B' the X-ray binary jets would be injecting a 
significant amount of power into the ISM. 



CONCLUSIONS AND CAVEATS 

In this paper we have demonstrated that the integrated ki- 
netic power from X-ray binary jets in our galaxy is sufficient 
for them to be a significant source of cosmic rays. Observa- 
tions have already clearly established that jets from powerful 
black hole transient XRBs can accelerate electrons to very 
high (TeV) energies (e.g. Corbel et al. 2002), but there re- 
mains little direct evidence for the acceleration of protons 
or atomic nuclei to relativistic energies (although heating 
of atomic nuclei to ~ 10 7 K at multiple sites in the jets 
of SS 433 has been demonstrated e.g. Marshall, Canizares 
& Schulz 2001 and references therein; Migliari, Fender & 
Mendez 2002). However, unequivocal evidence for acceler- 
ation of the non-leptonic component of cosmic rays in su- 
pernovae remnant shocks was also very thin on the ground 
until recently (Aharonian et al. 2004), and historically the 
argument was based upon energetics alone. 

Given that the underlying physics of particle accelera- 
tion in SNR shocks and XRB jets / jet-ISM interactions is 
likely to be similar, if one process can accelerate protons 
to very high energies then presumably so can the other. 
It should furthermore be noted that XRB jets almost cer- 
tainly spend a considerably longer phase moving with bulk 
relativistic velocities than do SNR, and the variation of effi- 
ciency of shock acceleration with bulk Lorentz factor is not 
well known. So it returns, at this stage, to simply a ques- 



tion of energy budgets. Within the Milky Way, we conclude 
that the integrated jet power is likely to be around 1% of 
the rate of injection of energy from supernovae, which is 
around 5-10% of the cosmic ray luminosity, comparable to 
the numbers estimated by Heinz & Sunyaev (2002). There 
exists a distinct possibility that these numbers are serious 
underestimates. In addition, the total number of black holes 
in our galaxy accreting at very low levels remains highly un- 
certain, and the contribution from such sources is difficult 
to quantify; this is an area which needs some consideration. 

The situation in the early-type galaxies, as illustrated 
with the example of NGC 4472, is therefore even more strik- 
ing. In this system, with a mass 3.5 times that of the Milky 
Way, the total supernovae rate is likely to be only half that of 
our galaxy, yet with a a population of LMXBs which scales 
as the mass. Even if the energisation of cosmic rays in our 
own galaxy is dominated by SNe, in such systems a signf- 
icant fraction of the cosmic ray energy input should come 
from X-ray binary jets. 

Furthermore, by assuming, as suggested by Gilfanov 
(2004), that the normalisation of the LMXB X-ray lumi- 
nosity function is proportional to a galaxy's stellar mass, we 
have most likely underestimated the importance of LMXBs 
in giant elliptical galaxies and especially in the central dom- 
inant (cD) galaxies of clusters; some scatter is found in the 
relation of Gilfanov (2004), and in fact, the residuals are cor- 
related with the globular cluster specific frequencies of the 
galaxies (White, Sarazin & Kulkarni 2002). Furthermore, it 
is known that globular cluster specific frequencies are gen- 
erally higher in the most massive galaxies (e.g. McLaughlin, 
Harris & Hanes 1993), and the metal rich globular clusters, 
which are more likely to have X-ray binaries (Kundu, Mac- 
carone & Zepf 2002) are even more prevelant in these sys- 
tems (e.g. Kundu & Whitmore 2001) - the fraction of the 
X-ray sources in globular clusters clearly increases as one 
moves from SO to giant elliptical to cD galaxies (Maccarone, 
Kundu & Zepf 2003). As a result, the most massive galax- 
ies are likely to have somewhat higher Llmxb / 'Lnir ratios 
than the Milky Way, making the estimates here conserva- 
tive. 

There are of course many caveats and unknowns asso- 
ciated with the entirety of the above discussion, as there are 
with the assumption of SNe-dominated energisation. These 
include the energetic contribution from gamma-ray bursts 
and 'hypernovae' (e.g. Iwamoto et al. 1998), although the 
rates for such events are likely to be orders of magnitude 
less than those of core collapse supernovae (Podsiadlowski 
et al. 2004). Furthermore, there exists a distinct possibility 
that there is a poorly-known population of radiatively-faint 
sources with powerful jets such as LS 5039 (Paredes et al. 
2000). 

To conclude, we argue that the input of kinetic energy 
into the ISM in our galaxy may have a significant contribu- 
tion from the jets of X-ray binaries, although is still likely to 
be dominated by supernovae. In early-type galaxies, which 
have no core collapse supernovae but a comparable number 
of low-mass X-ray binaries per unit stellar mass, the relative 
contribution will be much larger and should be considered 
as potentially rivalling supernovae. Given the high Lorentz 
factors of the most powerful X-ray binary jets, their effect 
may be most evident in its contribition to the production 
of cosmic rays. Finally, the energetic contribution of these 
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sources to the ISM will be in the form of a much more dis- 
tributed source than the rare, more energetic, supernovae, 
which may have implications for the distribution of cosmic 
rays throughout the galaxy, and be an important input for 
models of their propagation through the galaxy (e.g. Strong, 
Moskalenko & Reimer 2004). 
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